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Abstract 

We investigate the cosmological perturbations in generalized gravity, where the Ricci scalar and 
a scalar field are non- minimally coupled via an arbitrary function /(</?, R). In the Friedmann- 
Lemaitre-Robertson- Walker (FLRW) background, by studying the linear perturbation theory, we 
separate the scalar type perturbations into the curvature perturbation and the entropy pertur- 
bation, whose evolution equations are derived. Then we apply this framework to inflation. We 
consider the generalized slow-roll conditions and the quantization initial condition. Under these 
conditions, two special examples are studied analytically. One example is the case with no entropy 
perturbation. The other example is a model with the entropy perturbation large initially but 
decaying significantly after crossing the horizon. 
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I. INTRODUCTION 



In the past three decades, huge progress has been made on our understanding of the 
early universe, both theoretically and observationally. This is implemented by the inflation 
theory [J 0, |jj merging the general relativity and quantum field theory in an elegant way. 
On the one hand, inflation theory has naturally explained the initial condition of big bang 
cosmology. On the other hand, it also makes quantitative predictions which can be tested 
by precise observational data 0, 0, 0] . 

So far the prevail inflation models are based on the Einstein gravity coupled minimally 
to a scalar field (or more scalar fields) @, 0, fi M, 

I 111 . Il2l . 1 131 Whereas considerable 



investigations were also performed on models of modified gravity with higher derivative 
corrections or non-minimal coupling. Most of them can be classified into two categories: 



• f(R) models without a scalar field 15. 16l]; 

• F(<p)R scalar-tensor theory. 

Each of them has only one degree of freedom. This is clear if one takes a conformal trans- 
formation as done by 0,0,0- 



Our intention in this paper is to deal with a general type of model, namely f(<p,R) 
gravity, which unifies and generalizes the above relatively simpler models. The action of this 
model is of the form 1 



S = d 4 xy/^g 



l -f(^R)- l -g^d^d^-V{ V ) 



(1) 



Here the f((p,R) term contains a non-minimal coupling between the scalar field <p and the 
Ricci scalar R. While V(tp) is the potential of the scalar field. In principle V(<p) can be 
absorbed in f(<p, R), but we will keep them separate. 

In contrast with simpler models, the f(<p, R) model usually introduces another degree of 
freedom. Generally speaking, due to the new degree of freedom, there is a non-vanishing 
entropy perturbation in most models based on f((p, R) gravity. In the present work, we 
will distinguish the entropy perturbation from curvature perturbation, and then study their 
evolutions. 

To make our discussion self-contained and clear in notations, first of all, we collect some 
previously known results in section [Til Our general result is presented in section IIII[ where 
we extract the curvature perturbation and entropy perturbation as well as their evolution 
equations. Applying this formalism to inflation, we study the generalized slow-roll conditions 
and the quantized initial condition in section IIV1 In section [V] some typical examples 
are studied under the slow-roll approximation. One example is the case with no entropy 
perturbation, including the simpler models with one degree of freedom we mentioned above. 
The other example is to add a g(ip)R 2 correction to Einstein gravity. Specifically, we study 
the g{(f) = \\ip 2 , V((f) = \m 2 Lp 2 model in the limit M 2 /ip 2 Xm 2 (p 2 /M 2 <C 1. Ignoring 
the coupling of perturbations inside the Hubble horizon, we find the entropy perturbation is 
large at horizon-crossing but decays significantly outside the horizon. Initially the correlation 
between the curvature perturbation and the entropy perturbation has been neglected under 



1 Throughout this paper, we employ the reduced Planck mass M p = \/\J&nG and set H = c = 1. 
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our approximation, but at the end of inflation they become moderately anti-correlated. 
We summarize and refer to some open problems in section IVII For reference and as a 
support to our canonical quantized initial conditions, in appendix |A] we collect the relevant 
results of a two-field model which is conformally equivalent to the f(<p, R) model. In the 
generalized gravity, since it is difficult to draw a clear borderline between gravity and the 
matter, there is ambiguity in defining curvature perturbation and entropy perturbation, we 
present a more traditional (but less tractable) definition of these perturbations in appendix 
iBl Complementary to section IHH details for deriving the evolution equation of entropy 
perturbation are relegated to appendix 



II. REVIEW OF PREVIOUSLY KNOWN RESULTS 

The cosmological perturbations in f((p, R) gravity has been studied actively by Hwang 



and Noh in [20|, |21|, [22|, [23|, [24J etc. But most of the investigations mainly concentrated 
on theories with one degree of freedom, including the f(R) theory and the scalar-tensor 
theory. For generalized f(<p,R) theory, the complete evolution equations of the first order 
perturbations were obtained in [24]], where the background dynamics were also shown. In 
this section, to make our discussion self-contained and clear in notations, we write down 
these results following the notations of [7|, except for that the metric sig nature is taken to 



be (— , +, +, +). All of the results collected here can be found in [24|, [25|. Throughout this 



paper, we will focus on the scalar type perturbations, working in the longitudinal gauge 



The tensor type perturbation has been addressed in 23]. There was a discussion of scalar 
type perturbations in [22|, though, for generalized gravity, their full evolution equations were 
obtained in [24]. In this section, we summarize these results in a self-contained way, and at 
the same time, set up our convention of notations. In the subsequent sections, we will take 
further steps to get some new results. 

The background is described by the Friedmann-Lemaitre-Robertson- Walker (FLRW) 
metric 

ds 2 = ~dt 2 + a 2 (t)dx 2 = a 2 (r)(-dT 2 + dx 2 ), (2) 

where t is the comoving time and r is the conformal time, with respect to which the deriva- 
tives will be denoted by a dot overhead and a superscript prime respectively. Later on, for 
the sake of convenience, we will also use a superscript "•" to denote the derivative with 
respect to t. Then in terms of the Hubble parameter H = a/a, the Ricci scalar 2 can be ex- 
pressed as R = 6(2H 2 + H). Throughout this paper, we will only deal with the flat universe. 
The results for closed or open universe are expected to be similar. 

For succinctness let us introduce the notation F = -Jj^f((p,R). In this paper we will 
concentrate on the case with F > 0, but it is easy to extend our results to the case F < 0. 
The fluctuation of the scalar field <p will be denoted by Sip. It is well known that the 
perturbations of the metric can be decomposed into three types: the scalar type, the vector 
type and the tensor type. In the scalar-driven inflation, these three types are decoupled 
with each other if we only consider two-point correlation functions. Here we are interested 
in the linear order scalar type perturbations, so we can treat them exclusively, not worrying 



When defining the Ricci scalar R = g^R^, we take the convention of Ricci tensor R^ u = d\T* 
du^x + r^rj^ - l^ K r£ A with the affine connections I*, = ^{d^ + d u9llK - d K g„ v ). 
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about the vector type and tensor type perturbations. Considering scalar type perturbations 
only, the perturbed metric takes the form 



ds 2 = -(1 + 2(j))dt 2 + 2adiBdtdx l + a 2 [(l - 2^)<S y + 2d i d j £}dx i dx j . 



(3) 



Here 5ij is the Kronecker delta function. We will mainly work in the longitudinal gauge, 
that is, choosing B = £ = 0. If necessary, one can easily recover all of our results into the 
gauge-invariant form with the following dictionary 0, 24j |: 



(f) _> (f)(9i) =( f) + 



[a(B - £')]' 



$ _> =i/j- a H(B-£'), 
^ S(p {gi) = 5p + ap{B- £'), 
5F -> = + aF(£ - £ ')• 

Corresponding to action flTJ), the variation of (p and gives 



6iS 



d A x£> 



+ / d*x 



-9 



+d a Fg a ^{d p g^ + - d u g^) - d^Fg a ^{d a g Pu + - d u g a/3 ) 



-gd x F(2g al3 g XK - g aX g^ - g aK g^ 



where the total derivative term 



<p5ip) 



+d K 



l -^g-d x F(2g^g XK - g*g v " - g^ K g vX )5g^ 



Then the generalized Friedmann equations are simply 

\p 2 + V - + 3H 2 F + 3HF - 3HF = 0, 
if 2 + 2HF - HF + F = 0. 



(4) 



(5) 



(6) 



(7) 
(8) 



Equation (jSJ) proves to be very useful and we will use it frequently without mention. The 
non- minimal coupling f(<p,R) brings a new term to the equation of motion for the scalar 
field 

Cp + 3Hp - \f, v + V itp = 0. (9) 



From (jSJ), we get the Einstein equations of linear perturbations 24], |25j] in the longitudinal 



4 



gauge 



1 1 



3H(ifj + H<j>) - -V 2 if> = ~^P, (10) 

i + H*—^ (11) 
SF 

^ + 3#(#0 + ijj) + H<P + 2H<\> + ^V 2 (0 - V) = (13) 



p 



in which <5p, Sq and <5p are defined as 



M 2 



Sp = ~-£-[<p6<p + -(-/,„ + 2V; V )^ - 3#5F + (3# + 3# 2 + — V 2 )5F 

+ (3HF - 2 )0 + 3F(H(f> + if})), 
M 2 

= -^-di{ipStp + SF-H8F-F<f>), 
b 

M 2 1 .... 9 

Sp = -^[ipStp + ~(f j(p - 2V >ip )5<p + 5F + 2H5F + (-H - 3H 2 - -^V 2 )SF 

-F<j>- (<p 2 + 2F + 2HF)cf>-2F(Hcf> + ip)]. (14) 

Equations (I10H13P follow respectively from the G®, G®, G^(i ^ j), G\ components of Einstein 
equations. The equation of motion for Sip gives a redundant relation. 

Using equations ( fTTj) and ( fT2l) to cancel Sip and 8F respectively in ( fTUj) . one obtains 



F(4> + + (HF + 3F- —)(</> + i>) 

+ \(3HF + 3F) - ^(HF + 2F) - ^V 2 ]<p 

ip cr 

+ [(AHF + HF — F) — —(HF — F) — ^V 2 }ip = 0. (15) 

ip a 2 

Inserting (TTTT) and (|T2l into the G\ component equation (|T3|) . it will result in a relation 
automatically satisfied by the background equation 

Remembering that 8F = FrSR + F jip Sip and F = FrR + F }ip ip, we can write (|T2|) in 
another form 

F{<f> - if;) = 2F tR [34> + 3#0 + \2Hi) + 6(2# 2 + H)<f> + \s/ 2 ((f) - 2^)] 

(X 

-^[F{4> + if;) + (HF + 2F)0 + (FF - F)ip]. (16) 

In minimally coupled model / = M 2 R, equation (fTB"j) reduces to the familiar relation <p = if). 

Comments are needed here. We start with five equations (four perturbed Einstein equa- 
tions and one equation of motion for Sip) of three variables (Sip, <f>, if}), so this system of 
equations appears to be over-determined. However, as we have just mentioned, the equation 
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of motion for 5ip is redundant, which can be derived from perturbed Einstein equations. 
In addition, one of the remaining four equations turns out to be automatically satisfied by 
background equations. This is also understandable if one recalls the Bianchi identity. Con- 
sequently, there are three independent equations with three variable at last. What we have 
done was just using one of the equations to cancel 5<p, and arriving at two equations (TT51) . 
(I16p with two variables 0, ip. 



III. CURVATURE PERTURBATION AND ENTROPY PERTURBATION 



The perturbation equations (TTo^) and (TLT)j) were used in 24, 25j to study certain f(tp,R) 
models with only one degree of freedom, where the effects of entropy perturbation have been 
neglected all the while. However, in the most general case, the entropy perturbation may 
play an important role in the early stage and then decay or translate into the curvature 
perturbation, leaving some observable signatures in fluctuations of cosmic microwave back- 
ground (CMB) and dark matter. 3 Disregarding the entropy perturbation would obscure 
many interesting phenomenological predictions, such as the residual entropy perturbation 
and the large non-Gaussianity. From this point of view, for generalized f(<p, R) gravity, it 
is important to take the new degree of freedom into consideration and study the entropy 
perturbation seriously. 

In this section, we will rearrange the perturbation equations (IT5|) and ([TBI) in order to 
decompose the scalar type perturbations into curvature and entropy components and to 
get their evolution equations. In subsequent sections, applied to inflation, the evolution 
dynamics of perturbations will be studied under the slow-roll approximation. 

For our following study, it is essential to notice that equation (fl5l) can be recast in the 
form 



n 



In 



2Fip 2 + 3F 2 J 2Fip 2 + 3F 2 
with the curvature perturbation 



2HF + F . 2HF + F F ~. ± 

Ss + -zr-, — — — V 2 (0 + VO, (17) 



2Fip 2 + 3F 2 a 2 



1 , s 2HF + F 
2 Vr V> 2F<p 2 + 3F 2 



F(<f> + + 2 ( 2HF + F )( ( f ) + VO 



(18) 



The first term on the right hand side of (|17[) can be taken as the entropy perturbation, 4 or 
more exactly, the relative entropy perturbation 13l 



F{2HF + F) l~Y §s 
(p{2F(p 2 + 3F 2 ) V 2F S ' 



(19) 



A mechanism for entropy perturbation decay into curvature perturbation was realized in the curvaton 
scenario 



2G 



27j . In /(</?, R) gravity one may expect a similar story. It is remarkable that by curvaton 



2G 



27 



28 



29 



30] 



mechanism a large non-Gaussianty can be generated 

Here is an ambiguity in normalizing entropy perturbation. We define the entropy perturbation by (Q3 
but we will abuse it for Ss and 5s since they are proportional to S up to background quantities. Through 
relations ([TO)) and (|28p they are easy to be traded with each other. Our normalization of (|19p is chosen 
to ensure that V-r* = 'Ps* when perturbations cross the Hubble horizon, as given by equation J72|). 
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where 



5s 



• + ^) + i(2#F + F)(0 + ^) + 



2F<p 2 + 3F 2 fA 



2F 



(20) 



Strictly speaking, the second term on the right hand side of ([I 



also contributes to the 
Throughout our 



total entropy perturbation, but it is suppressed on super-horizon scale [13 
paper, we focus on the relative entropy perturbation, and refer it as entropy perturbation 
for simplicity The adiabatic (curvature) and isocurvature (entropy) perturbations were 



investigated in [3jJ, |32j, |33|, [34j for a two-field Lagrangian with a non-standard kinetic term. 



The Lagrangian discussed there is equivalent to the f(<p, R) theory here, up to a conformal 
transformation [l8|, 2l[. One can check that the curvature perturbation defined in ffTSl) is 
conformally equivalent to that appeared in [3l|, 32, 33|, 34j. See also appendix fA] As another 
check, in the minimally coupled limit / = M 2 R, the expression ffTSl) reduces to the familiar 

form n = (j)- H(4> + H6)/H. 



Along the line of [35 



we give an apparently more traditional definition of curvature 
perturbation H e ff an d entropy perturbation Ss e ff in appendix [B] At the first glance, the 
traditional definition seems more physical. But it depends heavily on an artificial separation 
of gravity and matter content. The choice given by (|18|) and (120]) is more convenient in 
calculation. Moreover, the curvature and entropy perturbations at the end of inflation are 
not necessarily the ones probed by astronomical observations 0,0], because they may evolve 
significantly after the exit of inflation, depending on the details of reheating. We will take 
a pragmatic attitude and prefer the convenient definition (fl8l) and (1201) here. 

One may still raise a question: why do we claim that TZ and S defined above correspond to 
curvature perturbation and entropy perturbation respectively? This doubt can be resolved 
by rewriting ([TBI) and (1201) into the following form: 



K 



5s 



2HF + F r 3HF 

ipSip + 



2Fp 2 + 3F 2 
5<p 5F 



2F(p 2 + 3F 2 



-5F, 



(21) 



When deriving these relations, we have used equations (jHJ), ( fill) and ( JT2l) . 

First of all, with ( T2T1) at hand, we can apply it to the special limits mentioned in section 

m 

• f(R) models without a scalar field; 

• F(<p)R scalar-tensor theory. 



These models are extensively studied in the past [20|, |2jJ, [22|, [23|, [2J]. For f(R) models, 



0, thus the entropy perturbation vanishes obviously. As for scalar-tensor theory, since 



F becomes a function exclusively dependent of <p, we have 

5_F_ _ dip 



(22) 



which guarantees 5s = 0. Now it becomes very clear that there is no entropy perturbation 
in these models. This is in agreement with the fact that there is only one degree of freedom 
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in these models. On the other hand, curvature perturbation 7Z reduces to 



, ip + QSip, for pure f(R) gravity without scalar field; 
' V + j5F, for / = F(<p)R, ^ 



These are exactly the ones having appeared in {24|. Our definition of curvature perturbation 
1Z naturally generalizes them in a unified form. 

Secondly, in favor of the small dictionary (jlj), one can prove that the definitions of 1Z and 
5s are gauge-invariant. In other words, the (B — £') terms cancel out automatically. 

As a further check, since our definition of curvature perturbation is gauge-invariant, by 
choosing a special (but not longitudinal) gauge 

(2HF + F)v?<V c) + (3HF - <p 2 )5F {c) = 0, (24) 

we get a neat relation 

K = iJ {c) . (25) 

In minimally coupled models, the same relation (125]) holds in comoving gauge or uniform 
density gauge 5ip^ = 0. Here the gauge condition (l2"4l generalizes the comoving gauge 
condition, so we can take it as a generalized comoving gauge. In fact, just as in minimal 
models, relation (I25p has a geometric interpretation. The spatial curvature to the first order 
of perturbations is given by 

(3)^=1^). (26) 
cr 

Therefore, in gauge ff24l) . the adiabatic perturbation TZ is proportional to the spatial curva- 
ture ^R. This is why we give it the name "curvature perturbation". 

According to our definition in this section, the entropy perturbation is proportional to 
$}£_ _ fj^jg j g na t ura ] if one remembers that a new scalar degree of freedom is related to 

<p F & 

F in the present theory. It gets dynamical when F becomes dynamical. This form is very 
similar to the entropy perturbation in models with two fields, say (fi, <£> 2 , corresponding to 
an entropy perturbation proportional to ^ — 

We would like to emphasize that the perturbations defined in appendix [B] make senses only 
if one views the f(tp, R) theory effectively as Einstein gravity with exotic matter contents. 
This viewpoint is not so reliable since the gravity has been modified in f(<p,R) model. 
Furthermore, what we are really interested in is the evolution of (f> and ip. As we will show 
below, the choice in this section is powerful to study their evolution. 

In appendix O, we demonstrate that the evolution of entropy perturbation obeys the 
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equation 
5s -- 



In 



2 (2Fy? 2 + 3F 2 



+ • 



F 2 



, • V 2 
-3H\Ss+ -^5s 



-ip 2 -3F 



. ip 2 3F 
{ T + ^F 



,i , • •> •-. f " 2 FFF, 



2F(2H 2 + H) + — 7— - 2 - 3F - 



+ 



x • 



/ ln V9 2 (2F V 9 2 + 3F 2 ) y 



3# 



2F 



+ 



2F0 2 + 3F 2 
2Fy? 2 



5s 



3F 



+ F In 



In terms of 



2F0 2 + 3F 2 



5s 



2 3 •• 
^ + 2 F - 



V 2 

— (0 + ^)- 



3F R 2F j Rtp 
Fip 2 \' 3F(p 



this equation can be reexpressed as 



=5s, 



ipjAFtp 2 + QF 2 



. V 2 

<5s + 3if <5s ^-5s 

cr 



1 

2 



' y? 2 (2Fy? 2 + 3F 2 

In : 

F 2 



3, y? 2 (2F0 2 + 3F 2 ) 

- In : 

2 F 2 



-3H 



In 



V 2 (2F0 2 + 3F S 



F 2 



'5s 



+ • 



~2\ •• 



- y? 2 - 3F 



. y? 2 3F 
1 F + 2^ 



2F(2# 2 + if) + -<p 2 -3F 



FFF, 



3F R 2F R <p 



+ 



x 



In 



y? 2 (2Fy? 2 + 3F 2 ) 



2F 



2F<p 2 + 3F 2 



F 2 



^5s 



3H 



+ 



™£ + F|ln 

3F 



2F0 2 + 3F 2 



^2 \ F J 



V 2 

— (0 + ^)- 



ipjAFip 2 + 6F 2 



a 



(27) 



(28) 



(29) 
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At the same time, according to (fITI) . the evolution equation of TZ is 



K 



In 



2(p{2EF + F) 



5s + 



2HF + F F 



2F0 2 + 3FV F X / 4 F02 + 6 F2 2F0 2 + 3F 2 a 



V 2 (0 + ^), 



(30) 



There are similar equations in non-standard two-field models [3ll . |32| . l33l . |34J |. In prin- 



ciple, equations (1291) a nd (pJJ) ma Y be obtained by a conformal transformation from the 



counterparts in [33|, [34J. However, having avoided the intricacies of conformal trans- 

formation at the perturbation level, our derivation in the f(tp, R) frame is strai ghtf orward 



Nonetheless, it is still interesting to compare the results here and those in |3ll . |32| . l33l . 134 



Such a comparison would confirm the conformal 



via the conformal transformation [18|, [21 
equivalence at the perturbation level. 

Until now we have not made any approximation. Equations (129]) and (|30|) determine the 
dynamics of entropy and curvature perturbations exactly. They are applicable in various 
cosmological stages and scenarios of FLRW universe. Given a concrete model, the classical 
evolution of perturbations can be numerically followed utilizing these equations. To work 
in the inflation scenario and generate an appropriate spectrum of density fluctuation, we 
should consider the generalized slow-roll conditions and the initial condition. This is a task 
of the next section. 



IV. SLOW-ROLL APPROXIMATION AND QUANTIZATION 

In this section we will study the perturbations under the generalized slow-roll approxi- 
mation. Firstly we study the classical evolution. At the end of this section, we will discuss 
the quantization of perturbations as an initial condition. 

From equation (Tl8l) . we know that the curvature perturbation TZ is fully determined by 
cf) + ip and its time derivative. So the dynamics of <fi + ip informs us the dynamics of TZ. 
Inserting (Tl8|) into (|T7|) . one obtains 



In 



1F(p l + 3F 2 
¥ 2 



5s + (0 + ip) + In 



aF 3 



2F0 2 + 3F 2 



(0 + ^) 



+-[ln(a 2 F)]- 



In 



(2HF + F) 2 
2Fcp 2 + 3F 2 



V 2 



(31) 



We observe that substituting this equation into (1271) will result in a fourth order differential 
equation of <fi + ip. But that equation would be rather difficult to solve directly. In this 
section, we rewrite equations (12 7p and (1311) into a tractable form. In the next section, we 
will get their analytical solutions in some special examples. 
Firstly, let us define two new variables: 



fI 



:(0 + u s = a5s. 



(32) 



'4F0 2 + 6F 2 

Similar to 4>+ip, the variable u-ji tells us the evolution of curvature perturbation TZ. Following 
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relation (}3~Tj) . it obeys the equation of motion 



In 



with the "mass squared" 



m 



n 



2F(p 2 + 3F 2 



+i[ln(a 2 F)]« 



In 



2F(p 2 + 3F 2 
J 2 



2FCp 2 + 3F 2 



(2HF + Ff 



u s + u^ — V 2 m-^ + m 2 z a 2 UTi = 0, 



(33) 




2F(p l + 3F 2 
^F 5 



In 



2Fip 2 + 3F 2 
F~ 3 



2Ftp 2 + 3F 2 
If we introduce the notation 

E 

then one can easily prove that 

E = 



2HF + F 

t 2 



2Fp + 3F 2 
2F% 



(34) 

(35) 
(36) 



During the inflation stage, the Hubble parameter and the effective energy density are al- 
most constants. In the following, we will study perturbation dynamics under the slow-roll 
approximation. For the convenience of calculation, our definition of slow-roll parameters is 



H 

61 -If 2 ' 


62 = 


H 

HH' 


* = H^ 




712 = H<p> 


Si 


F 

HF' 


52 ~ HE' 




HF' 


s 4 


E 
' HE' 


A ^ 
06 = 

HE 


(37) 


While the slow-roll conditions are given 


by 








M < i, 


h 


«1, 


\Si\ < 1. 


(38) 



Be careful with the notation and sign difference between the slow-roll parameters here and 
those in most literature. Under the slow-roll conditions, the background equations (jTJlHD are 
significantly simplified, 



V- l -f + 3H 2 F 

,u2 



0. 



0. 



3H^ - -f iV 



V ~ 

V,tp — 



0. 



(39) 
(40) 
(41) 



(p A + 2HF - HF 
1 

T 

Notice that equation (j4Tj) can be consistently derived from the leading order Friedmann 
equations (1391) and (14*01) . Neglecting higher order terms, we have 

■- 2 3F(p 2 -4Fip0 , 1 



F- r 

b 2 
<5 4 ~ 2?7i 



E 



2^ 



2F\ 
Vi - 



ei - -Si, 



5 3(5^3 - <5i?7i - 4771^2 
2 1 35 l -4r ]l 



(42) 
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Under the slow-roll approximation, the Hubble parameter is almost constant during in- 
flation. In terms of slow-roll parameters, the nondimensionalized "mass squared" for u-ji 
is 

2e 1 - m . (43) 



H 2 

Therefore, to the leading order, the evolution equation of u-n is 

u nk + k 2 u nk + m 2 n a 2 u nk + (3u Sk = (44) 

in the Fourier space. Here we have used the notation 

P ~ siga((p)aH^i - 2e x (45) 

with sign(yj) = ip/\ip\. In other words, the sign of (3 depends on the sign of tfi. 

At the same time, in terms of us, the evolution equation (1291) of entropy perturbation 
can be written to the leading order as 

u 'sk + k 2u sk + m 2 s a 2 u Sk + ak 2 u nk = 0, (46) 



in which 



ml 5 r F FF V , , 

a F 



aH HipFi 



3F \ 2F(p 2 + 3F 2 



sign(0)^ v /<y 1 -2ei. (48) 



In the above, we discussed the classical perturbations and a generalized slow-roll con- 
dition. The initial condition is governed by quantum theory of fluctuations. That is, we 
should expand action flTJ to the second order with respect to perturbations, taking the form 



5oS = / d 3 xdr 



\{d T v n ) 2 + \{d T v s ) 2 - l(d iVn ) 2 - \{d lVs f + Cjv n d T v s - V(v n ,v s ) 

(49) 

In this action the covariant variables v-r and vs are linear superpositions of perturbations. 
Then the canonical quantization leads to the initial condition 

v nk -> -^=e~ ikT e nk , v Sk -> -^=e~ ikT e Sk (50) 
V2k \J2k 



in the short wavelength limit k/(aH) — > oo. Here {e-jik^sk} is the orthogonal basis 

(e Q fc, epv) = S al3 5(k - k'), a,P = 1Z,S. (51) 
From the argument in appendix we believe the canonical variables in f(<p,R) gravity 



should be 

2F 



2F(p 2 + 3F 2 



<p6<p + + 3HF \ <f> + 10 

2HF + F 



2F 2F<p 2 + 3F 2 ^ 
2Fip 2 + 3F 2 2HF + F ' 

^ a [F6 V + F^{<p - ij)} = 6s = V6u s . (52) 



2Fcp 2 + 3F 2 cp \l 2F(p 2 + 3F 2 
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A detailed form of (l4~9l) and a closer analysis of canonical quantization will be presented 
elsewhere. 

When solving equations (jHJ) and (|46j) . the initial condition (I50I) for v-r and t>s sets the 
initial condition of ur and u$. In the short wavelength limit, we find v-Rk —> ^ u Rk- Thus we 
have 

u ™ - pij^ ,<(, '* )^ »■ ust ^7m e ~ ,tTist - m 

One can check that this solution meets (jHJ) and P61) well in the limit k/(aH) —>■ oo. 



V. ANALYTICAL EXAMPLES 

In equations (jSj) and (H6l) . it is clear that and Usfc are coupled, while the coupling 
strength is controlled by a and (3. In general the coupled equations are difficult to solve 
analytically, thus numerical approaches are needed. Nevertheless, under some circumstances, 
analytical results are available. Let us focus on the case 5s = in subsection IV A[ In 
subsection IV B| we discuss the g(ip)R 2 model with a non- vanishing entropy perturbation. 
Under certain approximation, we solve a special case analytically, resulting in nearly scale- 
invariant power spectra. 



A. No Relative Entropy Perturbation: 5s = 

This is essentially the case with only one degree of freedom. Therefore, when 5s = 0, 
equation PS"]) is expected to be satisfied automatically, while the dynamics of the survival 
degree of freedom is described by equation (jH|). Although we cannot give a general proof 
to this claim, we can check it in two classes of models mentioned previously: 

• f(R) models without a scalar field; 

• F(ip)R scalar-tensor theory. 

For f(R) models, this is quite clear because us = and a = 0. For F(ip)R scalar-tensor 
theory, it is less obvious. But looking at equation (|47|) . one may notice that we have assumed 
F,r 7^ when writing down (1461) and there is a singularity in m? s . Multiplying equation (14"6"|) 
with F t R to eliminate the illegal singularity in m|, we will find the resulting equation is 
well-defined and automatically satisfied. 

We can not exclude the possibility that the above claim does not hold in some cases, 
although so far we can not come up with such an example. In those cases, equation f|46|) 
gives u-ji = 0, which also satisfies (|44"|) . But then there will be no quantum fluctuations. So 
we are not interested in this possibility even if it exists. 

Hence let us assume equation PS"]) is satisfied more generally in the cases we are interested 
with 5s = 0. Then we are left with one differential equation of un, 

u nk + k 2 u nk + m\a 2 u nk = 0, (54) 

which is obtained by setting u$ = in equation f|44|) . 

In accordance with initial condition fl5^|) . we find the solution for (I5"4"|) is 

unk = -^je^-^V^tfH^i-kr) = C^zH^\z), with v 2 = \ - ^| . (55) 
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Here we introduced the notation z = —kr. This form of solution is consistent with the 
minimal inflation model For u ~ |, it gives a nearly scale-invariant power spectrum. 
The coefficient C is determined by quantization condition (1531) . Keep in mind that in the 
limit z — > oo, 

H M(z) - a /Ae^-^-i) oc e"^. (56) 

V 7TZ 

According to inflation theory, the CMB fluctuations are seeded by primordial quantum 
fluctuations inside the horizon in the early universe. When the primordial fluctuations 
crossed the Hubble horizon, they began the classical evolution phase. So the initial condition 
for classical evolution is dictated by quantization. Here we did not study the quantization 
condition carefully. It would be important to put this result on firmer ground by quantizing 
the second order action directly. 

If the scale factor is exponentially growing a ~ e Ht , then there is a relation 

aH = -- = -. (57) 

T Z 



Straightforward calculations give 



1 - 2v 



u" 



Rk 



I4~z 

Ce{yz+ l -^^H^\z). (58) 



Notice that a prime denotes the derivative with respect to conformal time r. 

In the long wavelength limit kr — > 0, for v ~ |, we obtain to the leading order of slow-roll 
parameters, 




2 




2HF + F\ I F f-kr 
2F(p 2 + 3F 2 



2Fip 2 + 3F 2 

As a result, the power spectrum of curvature perturbation is 



(59) 



k 3 . m l2 H A F . . 

VK = = + (60) 



and thus the spectral index is 



AH F 2Fy? 2 + 3F 2 ' , . 

n n - 1 = — + i — — — J— . 61 

H 2 HF H(2F<p 2 + 3F 2 ) 1 ; 



Recalling that in |24j Hwang and Noh gave the corresponding result 

|^| + — for pure f(R) gravity without scalar field; 

^ - 1 - S ML + J__ (2F^+3F*r for f _ p( )R ( 62 ) 
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Obviously the final results (IBTI) and fl62|) are perfectly matched. Notice that for the case of 
pure f(R) gravity we can recover fl62|) from (IBTj) by setting ip = 0. 

In the leading order, we also find n-ji — 1 = 1 — 2v. This indicates lZk/{HlZk) — for long 
wavelength perturbations. That is to say under the generalized slow-roll approximation, 
the curvature perturbation is almost conserved outside the horizon. On the one hand, this 
extends the previous result in [24J to the general case 5s = (and equation (l4"6"j) satisfied 
automatically). On the other hand, to get richer phenomena, we should take the entropy 
perturbation into account and solve evolution equations (jHJ) and fH6|) more generally. 

B. g(tp)R 2 Correction 

A relatively simple but non-trivial model with non-vanishing entropy perturbation is to 
consider the R 2 correction with a ^-dependent coefficient, 

f(v,R)=M 2 R + g(v)R 2 . (63) 

For this class of model, under the slow-roll approximation, equations (l39H4ip become 

V~3M 2 H 2 , if 2 + 2HF - HF ~ 0, 3H<p ~ 12g )ip H 4 - V v . (64) 

Then we get the following relations: 

9 9 8qV 

F = M 2 p +2gR^M 2 p +^-, 

p 

V tV <p ~ QM 2 p HH = 2HVe 1 , 

VJp + V w if 2 ~ 6M 2 (HH + H 2 ), 

3HV+6Hf + 3Ht*(*^-V, v ) <?• (65) 

^ P ;f P iW 

The slow-roll parameters (1371) can be expressed in terms of g and V and their derivatives 
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with respect to tp, 



H Ag„V„ M 2 p V 2 



H 2 Ml 2V 2 



^ 8g^V \2g„V v _ M 2 p V w M 2 p V 



2T/2 



Hip Ml Ml V 2V 2 ' 

F 16eiVXgV% £ 1 

<1i = ~ — — Oo = ~ f i 111 

1 #F " V, V (M* + 8gV) ' 2 1 2 

c ^ ,2e 1 V(gV), vv E 3 

63 = JFf = 7]1+ m ' ^hE* 2 *-**' 



Cl r]xHH VxV v V ^ 



x F 5 3<5i5 3 - 5i?7i - 4r]i?72 

* = SI = * " 2* 1 + • («) 

Subsequently, the coefficients in equations (jUJ and (JIB"]) are 

2 



m%, ~ # 2 (2ei - ?7i), (3 ~ sign(v2)a# V^i ~ 2e i> a - sign(0)-ai7A/(5i - 2e u 



4 

2 _ rr2 ' 



m 5 ~ if 



2 5l_3ei + 2^- 2+ gM 2 



(67) 



The above results are still too complicated to get some sense. Particularly, since there are 
so many slow-roll parameters, one may even worry about whether the slow-roll conditions 
( 158]) can be satisfied simultaneously. However, for the special case 5 

g{<p) = \\ip\ V{ip) = \m\\ (68) 
the slow-roll parameters take much simpler form as below: 

H 2AmV 2M p <P 8AmV 
Fi = ~ — r?i = — - — ~ — 

1 H 2 Ml <p 2 ' 11 Hip Ml ' 

F 4e x AmV E 1 

~~ ~HF ~ M p 4 + AmV' _ ~HE ~ 61 ~ 2 

h= H~F = 71l + 3Cl ' 6 ' = WE^ 2Vl ~ I* 1 ' 

H p 
£2 = —r-^Vu V2 = TT 7 - ^ + 3ei ' 



As a matter of fact, a model with R ip correction was discussed in 
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We have been concentrating on the case F > 0, so in this example we limit our attention to 
the case with A > 0. Obviously the slow-roll conditions (|38|) are satisfied when M 2 <C tp 2 <C 
M 2 /(Am 2 ). To meet this condition we should fine-tune the coupling constants to be very 
small Am 2 <C 1. This is the large field inflation. Although the value of field <p is larger than 
Planck mass, thanks to the small coupling constants, its energy density is still less than the 
Planck energy density 

In the above special case, if we further assume M 2 /(p 2 <C Xm 2 ip 2 /M 2 <C 1 during inflation, 
then we will find 

2AmV 

61 - — T/pi — ' m-h-^u 

p 

m\ ~ -2eiiJ 2 , (3 ~ sign(0)aif v / 2ei, 

2 

'3 

The model with coefficients ( F70T) is relatively simple. Let us discuss it in some detail. 

Again, due to the non-vanishing a and /3, the interaction terms form an obstacle to our 
analytical study. However, it is still interesting to make some rough estimates by neglecting 
these interactions before horizon-crossing/Hubble-exit. The problem is akin to the one we 
met in inflation of coupled multiple field. In fact, there is an excellent analysis of coupling 



| ~ (31ei - 2)H 2 , a ~ sign((p)^aHV2^. (70) 



effects in [37J. It turned out if the coupling terms are of slow- roll order, then they will give 
a correction of slow-roll order compared to the leading order contribution. In our case, the 
coefficients of coupling terms are of order 0(y/e[), so we expect their corrections to the 
power spectra are suppressed by 0(y/eT). 

Our scheme is taking limit M 2 /(p 2 <C Xm 2 ip 2 /M 2 <C 1, and disregarding the interaction 
terms related to a and /3 inside the horizon k/(aH) > 1. In accordance with initial condition 
(j53"|) . we get an analytical solution to equations (jUJ) and 



unu = -X^-toV^H£\-fr)ew,, with v\ = \ - 

u Sk = ^e^-^V^^H^i-kr)^, with v\ = - - . (71) 

From this solution, we find at the time of horizon-crossing, the power spectra of curvature 
and entropy perturbations are nearly scale- invariant, 



_ _ k 3 „ H 



k 3 , _ 10 H 



M 2 
p 



967r 2 A 2 m 2 ^- 



i 



27T 2 1 k * 47r 2 (f> 2 
while their cross-correlation is negligible, 6 



M 2 

(72) 



967T 2 A 2 m 2 (y9^' 



k 3 

V c * = 77^ (K k *,<Sk*) ^0. (73) 

Z7T Z 



6 This is because we have neglected the coupling terms in (|44[) and (j46|) before Hubble-crossing. The fact 
is, when taking coupling terms into consideration, we expect the cross-correlation is not negligible here, 
'Pc*/'Ps* ~ But to get the explicit value of it, one should either perform the higher order 

calculation or take a numerical method. 
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Therefore the spectral indices at the horizon-crossing are 



nn* - 1 = n s * - 1 = 4ei* - 2r) u = -Ae u . (74) 

The variables with an asterisk subscript take their values at the horizon-crossing time k = 
aH. Note we have chosen the normalization of S in (fT9"|) so that Vn* — Vs*- 

We would like to pause here and comment on the scale invariance of the entropy 
perturbation. If readers go through our calculation carefully, they would find the limit 
Mp/ip 2 <C \m 2 (p 2 /M 2 <C 1 is of key importance in making Vs* scale- invariant. In this limit, 
the Mp/(2gV) term is small, hence m 2 s /H 2 ~ —2 and subsequently vi ~ 3/2. If Mp/(2gV) 
is not small enough, we cannot get such a value of and then the entropy power spectrum 
will not be scale- invariant. Even in the limit we have taken, there are subtleties with the 
Mp/(2gV) term. If Mp/{2gV) ~ 0(ei), we should pick it up when calculating m|. How- 
ever, for simplicity of our following estimation, when writing down (j70p . we have assumed 
M}/(2gV)<£e 1 . 

If one considers interactions outside the horizon, it proves useful to describe the evolution 
of perturbations in terms of a general transfer matrix 

7?- \ / ! T ns \ ( JZ \ 



S) V T ss } \ S 

which gives the power spectra of scalar type perturbations at the end of inflation 

Vk = Vk* + T% s V s *, V s = Tl s Vs*, V c = T ns T ss V s *. (76) 

These are the power spectra probed by astronomical observation (1, @], unless they changed 
significantly from the end of inflation to the matter-radiation decoupling. 

Although the cross-correlation power spectrum is negligibly small when the perturbations 
cross the Hubble horizon, it may still be large at the end of inflation. This is is usually 
evaluated by the cross-correlation coefficient (3 [12J (here we use a tilde to distinguish it 
from the notation appeared in equation (1441) ) or the correlation angle A [14j : 

/3 — cos A — — ^=^=. (77, 



In terms of the entropy-to-curvature ratio 



" = & (78) 



it is simplified by the transfer relation (J76 



a ■ frr vr / r s* sign(T ss )T ns 
cosA = S1 gn(T 55 )T K5 ^ TT ^- = ^=_. (79) 

Here we have used the notation sign^ss) = Tss/\Tss\ and the fact r^* ~ 1. Following 
the line of 37J, we expect the coupling terms inside the horizon would introduce £?( A /e7) 
corrections to r^*. 

As a matter of fact, well after the Hubble-exit hj (aH) <C 1, the interaction term ak 2 u-jik in 
equation (j4"B"j) is negligible, thus the entropy perturbation evolves independently outside the 
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horizon. As usually done in literature, when dealing with the evolution equation of entropy 
perturbation, one can ignore the second order term (i.e. S k term) and terms proportional 
to k 2 . Recalling relations (fT9"|) . (1521) and 



aJF(2F<p 2 + 3F 2 ) 
u s = = S, (80) 

from equation (H6l) we can get the evolution equation of S to the leading order, 

S k = -j (j£ + 2 + e x - - 35 2 + ~5^\ S k = fi s HS k . (81) 

which is a first order differential equation because the second order term has been ignored. 
This equation, together with (|T7|) outside the horizon 



8 X V 3 V 2 
= A^ffS*, (82) 

dictates the transfer matrix (|75p . When evaluating m 2 s in equation f[8"Tj) . we should keep the 
slow-roll order quantities in (JI7]) because /xg is of the slow-roll order. Taking /x 5 and //^ as 
their average values between the horizon-exit and the end of inflation, we can quickly write 
down the solution for (IHTj) and (|82|) : 

S k = <S fc *exp Q" jj s Hdt^ = S k *e^ N *- N \ 

Tl k - Tlfc, = f linHSudt = [ ^S k *de-f^ dN = ^S k * [ e ^ N ^ - l] , (83) 
it, J Us l*s 

in which iV = \n[a en d/ a(t)] stands for the e-folding number from time t to the end of inflation. 
Using this solution one may check that S k /(HS k ) ~ 0(ei). So it was reasonable for us to 
ignore the S k term. 

With the above results and formulas, it is not hard to obtain the power spectra at the 
end of inflation: 

Vk^Vk+Vs.&^-W -l] 2 , 
Vs 

The spectral indices at the end of inflation are 

riT? — 1 = ns* — 1 n — -, 

H s + &[e»s(N*-N)-l} 2 

n s -l = n s *-l- 2/is, 

lis [2ew(".-") - 11 , 
n c - 1 = n 5 » - 1 eMJV .-A0 _ I ■ ( 85 ) 
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In model (J68]) with M^/ip 2 < Xm 2 (p 2 /M 2 < 1, we find 3if0 > 0, so sign(y?) = 1. Under 
the slow-roll approximation, we get 



/ 4ei 29 1 ^ 

If we estimate /ar and with their values at Hubble-exit, the final result will be very 
simple (taking N = at the end of inflation): 

V K ~ 1 I 12 f p -29e 1 »AT»/3 _ ^2 _^S_ _ -58 £1 »iV»/3 



V c 2 /3 



-29ei.JV»/3 / -29ei*iV*/3 



1) 



V s * 29 V e u 

46 17 

ra^ - 1 = -4ei*, n s -l = y^i*, ric - 1 = y e i*- ( 87 ) 

At the end of inflation, the curvature perturbation and the entropy perturbation are mod- 
erately anti-correlated, with the cross-correlation coefficient 



* =c " AH -*/iidriTi? (88) 



but the entropy-to-curvature ratio is small enough, 



Ts = 84K.. + 12 ■ (89) 

The evolution of power spectra of curvature and entropy perturbations and their correla- 
tion have been shown in figure [TJ We can see the power spectrum of curvature perturbation, 
denoted by the blue solid line, almost doubled from the horizon-crossing (at iV* — iV ~ 0) 
to the end of inflation (at iV* — N ~ 60). Although the entropy perturbation (denoted by 
the purple dashed line) was lager than the curvature perturbation at the Hubble-crossing, 
it was decaying exponentially with respect to — N in the super-horizon scale. The cor- 
relation between curvature perturbation and entropy perturbation is depicted by the brown 
dot-dashed line. Once crossing out the horizon, it kept a negative value, with its amplitude 
at first increasing and then decreasing. 

In figure [21 we illustrate the dependence of correlation coefficient (the top graph), the 
logarithm of entropy-to-curvature ratio (the middle graph) and tensor-to-scalar ratio (the 
bottom graph) on N* — N. It is clear that the entropy-to-curvature ratio dropped down 
quickly outside the horizon. The entropy perturbation and the curvature perturbation was 
totally uncorrelated (under the decoupled approximation) at the Hubble-exit but evolved to 
be moderately anti-correlated at the end of inflation. 

At the end of inflation, if one assumes n-ji — 1 ~ —0.04, then rs — 4 x 10 -6 . This is 
well below the upper bound put by the five-year Wilkinson Microwave Anisotropy Probe 
(WMAP) Of course, we cannot take these numbers seriously, because our above analyt- 
ical results are meaningful only for estimation. The coupling between entropy perturbation 
and curvature perturbation has been systematically ignored until horizon-crossing. In addi- 
tion we have estimated fi-ji and fi-ji with their Hubble-exit values. Remember that in writing 
down (J7D|) we have neglected the M^/{2gV) term. This may also have introduced some 
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Pr,s,c/Ps* 




Figure 1: The evolutions of power spectra with respect to e- folding number iV* — N 
after crossing the horizon. The curvature power spectrum V-r is signified by a blue 
solid line. The entropy power spectrum Vs is denoted by a purple dashed line. The 
cross-correlation power spectrum Vc is depicted by a brown dot-dashed line. All of the 
power spectra are normalized by Vs*-, the entropy power spectrum at horizon-crossing. 
The vertical black dotted line corresponds to iV* — N = 60. 

uncertainties if M*/{2gV) ~ 0(e\). To get a robust conclusion, one should take all these 
effects into account carefully, resorting to the numerical method. 

Another issue of our analysis is the special limit we have taken: M 2 /ip 2 <C \m 2 ip 2 /M 2 <C 
1. Taking such a limit is two-fold. On the one hand, it enables us to get a nearly scale- 
invariant entropy perturbation. On the other hand, we find 3Hip ~ m 2 if(Xm 2 ip' l —Mp)/Mp > 
in this limit, hence the scalar field was not rolling down but rolling up the potential V(ip). 
Also the Hubble parameter is forced to grow (e\ > 0). The same thing also happens in 
phantom" inflation (single-field inflation with a kinetic term of the wrong sign), see e.g. 



38j. In our model, as the inflaton ip grows, the slow- roll parameter ei tends to order unity 
and the inflation is expected to cease then. It is still unclear how to terminate inflation and 
trigger reheating in this model. We leave it as an open problem for future investigation. A 
more interesting problem is to get a more realistic, slow-roll model. 



According to the result in 24], to the leading order, the power spectrum of tensor type 



perturbation in f((p,R) gravity is 



V T ^ ■^ = V s *(25 1 -4e 1 ). (90) 
It is conserved after Hubble-exit, with a spectral index 

nT ^^-WF =2ei ~ 51 - (91) 
For the specific model (IBS]) , in the limit M 2 /ip 2 < \m 2 ip 2 /M 2 < 1, they are 

^L~4eu, n T ~-2e u . (92) 
rs* 

In our case, the tensor spectrum is red tilted, which is different from "phantom" inflation. 
One would have noticed that the entropy-to-scalar ratio is small at the end of inflation, 

^ 33644 , QQ2 (93) 

V n 841ei* + 12 1 ; 
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Figure 2: The evolutions of correlation coefficient cos A (the top graph), entropy-to- 
curvature ratio r$ (its logarithm, the middle graph) and tensor-to-scalar ratio (the 
bottom graph) with respect to e-folding number N* — N after crossing the horizon. 
The vertical black dotted lines correspond to iV* — N = 60. The horizontal black dotted 
line corresponds to cos A = —1, that is, the totally anti-correlated situation. 



Its evolution curve outside the horizon is plotted in the lower graph of figure El 



VI. SUMMARY 



We investigated the cosmological perturbations in generalized gravity, where the Ricci 
scalar and a scalar field are coupled non-minimally by an arbitrary function f((p,R), with 
the Einstein gravity as a special limit. This general form unifies the usual modified gravity 
39 and scalar-tensor theory, but often introduces an additional degree of freedom. In the 



FLRW background, by studying the first order perturbation theory, we decomposed the 
scalar type perturbations into the curvature perturbation and the entropy perturbation, 
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whose evolution equations are obtained. The effects of entropy perturbation in this class of 
model were seldom studied in the past. 

Then we applied this framework to inflation theory. The slow-roll conditions and the 
quantized initial condition are discussed. The quantization of second order action will appear 
in a future work, which will put our discussion at the end of section IIVI on a firmer ground. 

Analytically we studied two special examples: the 5s = models in subsection IV Al and 
the g(ip)R 2 model in subsection IV Bl In the former example we unified the previous known 
results in a unique form and extended them to the general case. The latter example is 
more interesting. In the model g((p) = |A<£> 2 , V((p) = ^m 2 {p 2 , we paid attention on the 
limit M 2 /(p 2 <C Xm 2 ip 2 /M 2 <C 1. This limit helps us to get a nearly scale-invariant entropy 
perturbation. Although the entropy perturbation was large at the Hubble-exit, it decreased 
significantly outside the horizon. At the end of inflation, the "entropy-to-curvature ratio" 
we defined by ( 1751) was of order 10~ 6 , well below the constraint of five-year WMAP data. 

In spite of the above results, there are still some important problems unsolved. First, 
the quantization should be performed to confirm the normalization in ( 1531) directly. Second, 
to get robust conclusions and richer phenomena, maybe one has to employ a numerical 
method to solve evolution equations (H4|) and (l4"6l) more generally. Third, the formalism we 
developed in section [Till can be applied to other cosmological stages and scenarios, such as 



the preheating stage and the curvaton scenario, etc. In |40|, |41|, |42j, |43| , different schemes of 
modified gravity were investigated mainly concerning their effects on the late time evolution 
of our universe, the possible effect of /(</?, R) gravity on late universe is still waiting for us 
to study. 
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Appendix A: TWO-FIELD MODEL AND CONFORMAL TRANSFORMATION 



For comparison, let us collect here the well-known results for two-field inflation with a 
non-standard kinetic term, which is conformally equivalent to the f(<p, R) generalized gravity 
17 . 18 . TgL 24]. Usually the /(</?, R) formulism is called Jordan frame while the non-standard 



two filed formulism is called Einstein fame. In this appendix, we use the super/subscripts "J" 
and "E" to distinguish different frames. The action two-field inflation with a non-standard 
kinetic term in Einstein fame is 



S 



d x^-g E 



1 



16nG 



Re ~ \gfd E aX df X - l -e 2 ^gfd E ipd E 



(Al) 



The reference [31] has an original and clear discussion on the curvature and entropy per- 
turbations in this model. The readers may get details there. In a further study [13j, it was 
clarified that the total entropy perturbation includes not only the relative entropy pertur- 
bation, but also an entropy perturbation proportional to k 2 /a 2 H 2 . Although it will give a 
correction at the scale k ~ aH, its contribution to entropy perturbation is suppressed in 
long wavelength, so we do not consider its effects in our investigation. 
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In longitudinal gauge, the canonical variable in this frame is given by 

v x = a E \ 5x + ~jf~^J , v v = a Ee b {Sip + -jj—^Ej • (A2) 

They are related to the canonical variables for curvature and entropy perturbations 

dfx e b dfip 



-K 



V(d t E x) 2 + e 2b (dfpJ" X VWW+^WW 2 * 



dfx e b dfy 
a E ^WW+^WW s (A3) 



Making use of the conformal transformation 



and the identification 



C = (A4) 



wrf^w (A5) 



one can prove the following relations jl8|, [2l|, \24 



dp, — (if, dtp — eft /, Hp, — — ~ 5 ■ — ~. 

E M p J E M p E 2F§ 

1, , - fs 



i>E = -((pj + iJj), S X = \j-M p {ijj-<Pj). (A6) 



If we take 



Kx) = -X/(VGM P ), Vfax) = ^[RjF(ip,Rj) - f(ip,Rj)}, (A7) 

then action ©can be perfectly reproduced from action (1A11) . For details of derivation, 
see reference [18]. In these frames, the conformal time are coincident cLte = drj, so the 
canonical quantized variables (1A3j) are exactly equivalent to (|52|) by the above conformal 
transformation. This conformal equivalence is powerful. Given a specific function f(ip,Rj), 
one can know the detailed form of V(cp,x) using F = j^jf((p,Rj) and (1A5I) . and then do 
our job in the more familiar Einstein frame. 7 



7 Thank the referee for an emphasis on this point. 
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Appendix B: A TRADITIONAL DEFINITION OF CURVATURE AND EN- 
TROPY PERTURBATIONS 



In all of our discussion, we take the curvature perturbation and entropy perturbation as 
defined in ( 118)) and f[2"U|) . But such a definition is different from the traditional one. If we 
regard the f{(p, R) theory as the Einstein gravity with exotic matter contents induced by the 
non- minimal coupling, and following the spirit of 35j, then we will arrive at a traditional 
form of curvature and entropy perturbations. Let us elaborate a little on this point. One 
should keep in mind that this point of view is different from that in appendix |A] Although 
we also use the term "Einstein gravity" here, it does not mean the Einstein frame. 

Formally, Friedmann equations (J7|) and (jSJ) can be rewritten as 



1 



3M2 



P, 



H 



1 



2M p 2 



(P + P), 



(Bl) 



with 



P 
P 



Ml 1 „ 
-f^ 2 + V + 
Ml ,\ 



p <-u? 
2^ 



V 



RF-f 
2 

RF-f 



-3HF), 

+ F + 2HF). 



(B2) 



Here in the effective energy density and pressure we have reckoned the contribution of the 
non-minimal coupling. Our final result seriously depends on this trick. The perturbations 
obey equations (U0H13p . Then the comoving curvature perturbation is given by 



H 



H 



The curvature perturbation on the uniform density hyper-surface is well defined, 

1 V 2 



(B3) 



H H / 

c e // = -^--^=-v>+^+#0, 

At the same time, the entropy perturbation 5s e ff is defined by 

T5s eff = 8p- c 2 s Sp, 

2 dp p 3HH + H 



1>. 



(B4) 



c: 



dp p 



It would be interesting to notice that 



Ceff + Keff 



2M p 2 



3HH 



(B5) 



3(p +p) a 2 



(B6) 



Finally, one can quickly prove 



H 
H 



H(j>) + 2H<p 



(B7) 
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and 



(B8) 



Now we see that the curvature perturbation 72. e // and entropy perturbation Ss e ff are 
related in the traditional manner. However, the full expression of <5s e // is rather complicated, 
accordingly its evolution equation is even more difficult to get. Therefore, the perturbations 
presented in this appendix are not convenient in studying the generalized gravity. Although 
the definition here appears very natural if one takes f(<p, R) gravity as an effective Einstein 
theory, a more convenient choice should be ([TBI) and (I20I) . 



Appendix C: EVOLUTION OF ENTROPY PERTURBATION 

From (fTgj) and (EDI), we get 



F(j> + VO = - \{2HF + F){cj> + V) + 2F ^ 2 + 3F ' 
2 V Mr ry 2HF + F 



<f>-ifj 



2F 



^5s 



2F 



2Fip 2 + 3F 2 2HF + F [2 
Now equations (j!5j) and (1161) can be rewritten in the form 



F(4> + ^) 



2F(p 



- HF -3F) (<p + ijj) 



+ 



+ 



^(2HF + F) — (2HF + FY 

3Fip 



+ 



(CI) 
(C2) 



(C3) 



F{4>-$) 



( AHF - 3F + ^ - (i + fy- 3HF(<j) - ^) 



+ 



+ 



2F(2H 2 + H)- (2HF + F)' + (2HF + F) 



2F(2# 2 + H) 



3F(p 



ip 2 -3F 



3F 



R 



V FF V \ 
FFF 



2F R( p 



(<f> + i>) 



2F, 

Iki 2 



+ ^) + -(0-V^). 



(C4) 
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Differentiating (!20j) with respect to time and making use of (1C3|) . one obtains 

'2F(p 3 



8s 



^)(0 + ^) + 2F ^l 3F2 (0-^) 



2F 



+ 



+ 



%{2HF + F) — -(2HF + F)' 

if 2 



— i* 



+ 



3 ~ 

2" 



(0 + ^) 



+-V 2 (0 + ^), 



which gives 



(0-^) 



2F0 2 + 3F 2 
2F 

/3 • 2Fif 

Ss+ [ -F - 

V2 



+ V0 + 



-(2#F + FY - ^(2HF + F) 

2 (p 



+ 



^)--v 2 (0 + vO- 



Taking the derivative of (1C5|) once more, we have 
6s 



f2Fp _3 p 



•■ 2F0 2 + 3F 2 ... 7 . 
2b 



+ 



+ 



+ 



2F(p 




+ —(2HF + F) — (HF + 2FY 



(0 + ^) 



^2F0 2 V | 3Fy3 



(0-^) 



+ 



-(2HF + F) 
F^" 



-{2HF + F)" + 



3Fy3 . 2 3 

cp i< 

y Y 2 

+^V 2 (0 + $ + (F - 2HF)^-(<P + V). 
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Substituting P2) and (lC4l) . flCHl) . flCil and (|C2]) into (E2]) step by step to eliminate ij> + <ij) 
and <p — ip, 4> — ip, 4> + 4> an d — V* respectively, one will arrive at a rather scattering equation: 



5s 



In 



<p\2F<p 2 + 3F 2 
F 2 



V 2 

3H )■ 5s + —5s 
a z 



+ < 



^(2HF + F) 





1 

2 1 



(21 


3F\ 




~2FJ 



^(2HF + F) - (2HF + FY 



+ 



+ • 



+ 



+ 



x • 



+ • 



+ 



+ 



In 



2F(2H 2 + H) - (2HF + /'•')* -i (2HF-\ F)(%- 

<p 2 (2F<p 2 + 3F 2 



FF, 



F 2 



3H 



-(2HF + F)* - %(2HF + F) 

2 Lf 



(f 6F R <p 



2F<f\' (p 



2(p 3F\ (2FCp 



+ ^(2HF + F) — (HF + 2F)' + - _ 



2F 



HF -3F 



<P 3F\ / „ ■ 2F(p F 2 F V \ 

+ — AHF - 3F + — £ 

F 2F \ <p 3F R <p) 



In 



2 (2F0 2 + 3F 2 



F 2 
2HF + F 



3H 
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2F 

Ffi 
~F 

<f 3F_ 
F + 2F 



3j, 2Fy3 
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2F 



F 2F / 2HF + F 



^2 U 2F 



3F/3 
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FFF, 
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3H 
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2 if 2 3F 



-F 



In 



2 (2Fy9 2 + 3F 2 

T 2 



] v 2 

3HF)— + V . 
I cr 



(CS 



This equation looks terribly lengthy. However, repeatedly employing the relation (jSJ) or 
namely 

(2HF + F)' = 3HF- tp 2 , (C9) 
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after careful calculation, we find the coefficients of the 2 ^p + p + 4>) — 7i] and remaining 
(<ft + ip) terms are exactly vanished, resulting in a much simpler form (j2"Tj) . 
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